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ABSTRACT 

O 
O 

The Gil. 11-0. 12 Infrared-Dark Cloud has a filamentary appearance, both in 

S— I ' 

absorption against the diffuse 8/zm Galactic background, and in emission from 
cold dust at 850/im. Detailed comparison of the dust properties at these two 
wavelengths reveals that standard models for the diffuse interstellar dust in the 
Galaxy are not consistent with the observations. The ratio of absorption coeffi- 
cients within the cloud is k 8 /k 85 o < 1010, which is well below that expected for 
the diffuse ISM where Ks/ K 850 ~ 1700. This may be due to the formation of ice 
mantles on the dust and grain coagulation, both of which are expected within 
■ dense regions of molecular clouds. The 850/zm emission probes the underlying 

O 

m 



radial structure of the filament. The profile is well represented by a marginally 
resolved central region and a steeply falling envelope, with S(r) oc r~ a , where 
r£3 • a > 3, indicating that Gil. 11-0. 12 is the first observed filament with a profile 

similar to that of a non-magnetic isothermal cylinder. 

O 

Subject headings: ISM:clouds - ISM: dust, extinction - ISM:evolution - ISM:structure 



ISMdndividual objects:Gl 1.11-0.12 



1. Introduction and Observations 

Infrared-dark clouds (IRDCs) have been identified by their substantial mid-infrared (8- 
25/xm) extinction in Midcourse Space Experiment (MSX) images (Egan et al. 1998). Egan et 
al. concluded from the high mid-infrared opacities of the IRDCs that they possess hundreds 
of magnitudes of visual extinction and contain large column densities of cold dust. They 
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derived characteristic physical parameters for the gas inside the IRDCs: I m 0.5 - 15 pc, 
T K ~ 10 - 20 K, n(H 2 ) > 10 5 cm -3 , and H 2 column densities ranging up to 10 23 cm -2 (Carey 
et al. 1998). 

A number of IRDCs have been mapped with SCUBA on the JCMT (Carey et al. 2000; 
Redman et al. 2000, 2003) and are found to have the mid-infrared absorption closely as- 
sociated with filamentary and/or floculent clouds in emission at 850 and 450/im. Bright, 
compact sources with masses in the range ~ 10 to 10 3 M and infall/outflow signatures are 
located along the filaments and appear to be in a variety of the early stages of star formation. 
The larger IRDCs appear to be similar in size to the Orion molecular cloud but different in 
many important aspects. They are more heavily extincted, colder, more quiescent, and at 
earlier stages of evolution. 

For the work reported in this paper, we investigate the IRDC Gil. 11-0. 12, at a kinematic 
distance of 3.6 kpc (Carey et al. 1998), because it contains several long, filamentary segments 
seen against a relatively uniform galactic background at 8/im (Fig. 1). Carey et al. (2000) 
first reported JCMT observations of this source. Star formation appears to be more active in 
the southern part of the cloud than the northern half, with bright, compact SCUBA sources 
(labelled PI, P6 and P7 in Fig. 2). The PI source is driving a molecular outflow and is 
associated with a weak point source of 8/zm emission. There is no evidence of high-mass 
star formation throughout most of the cloud. The catalog of Becker et al. (1994) lists two 
UCHII regions in the vicinity of Gil. 11-0. 12, but they are far enough away that they are 
unlikely to be affecting the IRDC properties even if they are at the same distance. 

JCMT 3 SCUBA observations were obtained at 850/im under fair weather conditions 
(t"850 = 0.24 - 0.30) in April 1999 with an angular resolution of 14.5". "Scan mapping" 
was employed to produce Nyquist-sampled images. The data were taken with the standard 
observing set-up using chop throws of 20", 30" and 65" at position angles of 0° and 90°. The 
data were processed using the JCMT data reduction pipeline (ORAC) to flatfield, extinction 
correct and calibrate the chop data (Economou et al. 1999; Jenness & Economou 1999). The 
calibration accuracy is good to ±15% (Jenness et al. 2002). The optimal map reconstruction 
(see Fig. 2) was then performed from the individual chop measurements, using a matrix 
inversion technique outlined by Johnstone et al. (2000). Slowly varying spatial features 
with wavelengths much longer than the maximum chop, 65", are not measured reliably; 
thus, large, low-amplitude structures in the reconstructed map tend to be artifacts. Further 
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descriptions of SCUBA and its observing modes can be found in Holland et al. (1999) and 
Jenness et al. (2000). 

Heterodyne observations of 13 CO (2 - 1) and C 18 (2 - 1) were taken near the P2 region 
of Gil. 11-0. 12 with receiver RxA3 on the JCMT during several observing sessions in July 
2000, April 2001, and June 2001. The velocity dispersion within the cloud can be estimated 
by fitting single component Gaussians to the line profiles. Typical linewidths obtained have 
a = 1.2 km s^ 1 for the 13 CO line and a = 0.9 km s -1 for the C 18 line. The temperature of 
the gas in the segment of the filament labelled P2 (Fig. 2) has been estimated from H2CO 
observations by Carey et al. (1998) to be less than 20 K. 

2. Dust Properties: Extinction versus Emission 

The extinction of background emission at 8/im is due to the presence of dust within 
the Gil. 11-0. 12 filament. This same dust radiates at submillimeter wavelengths and is 
observable as emission in the 850/im image. It is expected that there should be a strong 
correlation between the extinction and emission maps, and this is indeed seen in Figs. 1 and 
2. The correlation is complicated by several factors: (1) the diffuse background emission 
at 8/im, which the molecular cloud extinguishes, is not entirely smooth; (2) star formation 
occurring at PI produces 8/im emission; (3) the zero submillimeter emission baseline is 
poorly defined; (4) the measured emission and absorption properties are beam-convolved; 
(5) the dust temperature and emissivity may change with location in the cloud because the 
degree of dust coagulation and freeze-out of molecules onto grain mantles depend on the 
local physical conditions, such as density and temperature. 

Statistically relevant results can be obtained by analyzing a large sample of measure- 
ments, despite the above caveats. Fig. 3 shows the 8/im emission versus the 850/im emission 
for each location across the filament, along with the median 8/xm flux as a function of 850/im 
emission. The trend in both the data and the median curve is clear: the 8/im flux decreases 
due to extinction as the 850/im flux increases. The highest 850/im emission points do not 
follow this trend because of the presence of an energetic protostellar source (PI) within the 
deepest part of the cloud that is altering the temperature of the dust and producing weak 
8/im emission. 

Given the clear trends in Fig. 3, a simple dust model may be fit to the data. At each 
location x across the filament, the measured 8/im flux density /s( x ) depends on an unvarying 
foreground /s(fg) and background /s(bg) emission, and the optical depth through the cloud 
r 8 (x) = /t 8 S(x), where k x is the absorption coefficient of the gas plus dust at wavelength 
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X(/im) and S(x) is the column density of material within the intervening cloud. Explicitly, 

/ 8 (x) = / 8 (bg) exp [-« 8 E(x)] + / 8 (fg). (1) 

The column density of the cloud is related to the the observed flux density /s5o(x) (measured 
in Jy/beam), for the optically thin submillimeter emission, by the equation 

/s5o(x) = B 850 (T d ) k 850 £(x) tt, (2) 

where B\(T d ) is the Planck function at A and VL is the solid angle subtended by the telescope 
beam. Thus, the optical depth at 8/zm, which is given by r 8 (x) = k / 8 5o(x), is linearly related 
to the observed flux at 850/im, with k defined by k 8 /k 8 5o = k B SFM (T d ) fl. 

The correlation between the flux at 8/im and 850/xm in Fig. 3 allows for a determination 
of k (and / 8 (bg), / 8 (fg)). The overlaid dashed line in Fig. 3 is based on the best-fit model, 
with k = 3.75, / 8 (bg) = 2.3 x 1(T 6 Win" 2 sr" 1 , and / 8 (fg) = 2.5 x 10~ 6 Win" 2 sr" 1 . Excellent 
fits are found for a small range around k = 3.75 ± 0.25, implying 

k 8 _ 1350 
^ ~ exp(17K/T d ) - 1 

Observations of molecular lines provide a range of gas temperatures, T K ~ 10 - 20 K, with 
a best-fit value of T K = 15 K. Theoretical models of the dust temperature toward highly 
extincted cloud centers (Zucconi et al. 2001; Evans et al. 2001; Kramer et al. 2003) yield 
lower limits of 5 - 8 K. The corresponding limits on the opacity ratio are 1010 if T d = 20 K 
and 50 if T d = 5 K. 

For diffuse interstellar dust the expected opacity ratio is ~ 1700 (Li & Draine 2001); 
however, ice mantles form on the surfaces of dust grains in dense molecular clouds and 
the grains themselves coagulate. This increases the submillimeter opacity and obscures the 
silicate features that account for most of the opacity at 8/xm (Ossenkopf & Henning 1994). 
Ossenkopf and Henning found that, after 10 5 yrs of coagulation, dust with thin ice mantles 
entrained in cold (Tk < 20 K) gas with number density n = 10 6 cm -3 produced an opacity 
ratio between 8/im and 850;um of ~ 500, consistent with these observations if the dust 
temperature is T d ~ 13 K. Kramer et al. (2003) found a similar increase in the importance 
of submillimeter emission compared to (near) infrared extinction in the filamentary cloud IC 
5148. These result are consistent with evidence for dust growth in cold, dense protoplanetary 
disks where observations of the slope of the spectral energy distribution are used (Li & Lunine 
2003) 

All reasonable dust models require a foreground 8/zm emission of / 8 (fg) ~ /s(bg) which 
implies that half of the diffuse mid-infrared emission is produced in front of the filament. 
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This is intriguing given that Gil. 11-0. 12 is only 3.6 kpc away and on the near side of the 
Galactic Center. The diffuse 8/xm emission, while relatively uniform across the IRDC, must 
be significantly patchy on Galactic scales. 



3. The Filamentary Structure of Gil. 11-0. 12 

The simplest equilibrium model for a filament is that of an isothermal cylinder (Ostriker 
1964) where the density profile as a function of cylindrical radius r is given by 

P ~ f , . 9 /o 2 \2 ' (4) 

(1 + r 2 /8rg) 

where p c is the central density and r is the core radius defined by a 2 = 4nGp c r 2 . The 
velocity dispersion a should be derived from the temperature of the gas for a true isothermal 
cylinder. However, in most astrophysical cases the measured velocity dispersion within the 
molecular cloud includes a significant "non-thermal" (turbulent) component. 



The surface density profile for an edge on isothermal cylinder is: 

S 



<l +r 2/g r 2 



2N3/2' 



(5) 



where S = 2 1 / 2 7r p c r . The surface density profile becomes S « r 3 , at large projected 
distances from the center of the cylinder. 

Most studied filamentary clouds are not well-fit by this profile. L97 in Cygnus (Alves 
et al. 1998), the Integral Shaped Filament (ISF) in Orion A (Johnstone & Bally 1999), 
IC5146 (Lada et al. 1999), and a quiescent filament in Taurus (Stepnik et al. 2003) are all 
better fit by shallow density profiles with p oc r~ 2 , which are consistent with surface density 
profiles of S oc r~ l . Such profiles are well-fit by the Fiege & Pudritz (2000) model, which 
features a helical magnetic field. The r~ 2 density profile is obtained in FP models with a 
dominant toroidal field component, but steeper profiles are possible when the poloidal field 
dominates. Submillimeter polarization measurements trace the geometry of the magnetic 
field and provide additional evidence in support of the FP helical field model, which has been 
shown to agree with polarization maps of the ISF (Matthews et al. 2001) and NGC2024 
(Matthews et al. 2002). 

The radial structure of two well-defined regions of the Gil. 11-0. 12 filament, a northern 
segment including P3 and P4, and a southern segment including PI and P2, have been 
analyzed. The average intensity distributions orthogonal to the spine of the filament were 
determined using the same technique as applied to the ISF (Johnstone & Bally 1999), and 
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are presented in Fig. 4. The Gil. 11-0. 12 segments appear similar to each other, with a 
marginally resolved central region connecting to an outer region characterized by a density 
profile steeper than r~ 3 . This is in contrast to results for the ISF where the central region 
is unresolved, and whose surface density only falls off as r~ L0±a25 at large radii r S> r 
(Johnstone & Bally 1999). 

The emission from each section of Gil. 11-0.12 is reasonably fit by an isothermal profile 
(dashed lines in Fig. 4) allowing for an estimation of the core radius r . The core radius 
varies from r = 0.12 ± 0.02 pc in the northern segment to r = 0.10 ± 0.02 pc in the 
southern segment (assuming a 3.6 kpc distance to Gil. 11-0.12). The measured peak flux is 
fo = 0.20 ± 0.03 Jy/beam and f = 0.43 ± 0.06 Jy/beam, respectively, where the uncertainty 
in the peak flux is primarily due to uncertainties in the 850/im calibration. The observed 
dust emission can be converted to a central surface density S , despite uncertainty in both 
the dust temperature and emissivity, using 

S o = 0.13 f—^ ) f nn9 K85 ° 2 T ) 1 [exp(17K/T,)-l] gem" 2 . (6) 

\1 Jy/beam/ \0.02curg i y 

Assuming a typical dust temperature of 13 K, a gas-to-dust ratio of 100, and a dust emissivity 
K 850 — 0.02 cm 2 g -1 (Ossenkopf & Henning 1994), the central surface density is 0.07 ± 
0.01 g cm" 2 in the northern segment and 0.15 ± 0.02 gem -2 in the southern segment. If the 
dust temperature were much lower, the column density would be significantly larger. 

The central number density of H2 atoms can also be derived from the profile: 

n c = 2.6 x 10 4 ( £ ) ( ^-Y 1 ( ^—r) 1 [exp(17K/T d ) - 1] cm- 3 . (7) 

V 1 Jy/beam/ 1 VO.lpcy V°- 02cm S / 

For the canonical dust properties, n c = 1.1 ± 0.3 x 10 4 cm -3 in the northern segment, and 
n c = 3.3 ± 0.7 x 10 4 cur 3 in the southern segment. The measured C 18 linewidth observed 
for Gil. 11-0. 12 near the location of P2 is a = 0.9 km/s, allowing for an independent measure 
of the central density, n c = 3.4 ± 0.6 x 10 4 cm~ 3 , via the core radius definition. The agree- 
ment between these two derivations of central density provides support for the isothermal 
model. However, more complex filamentary models, including poloidal and toroidal magnetic 
components, may produce as good or better fits to the data (Fiege et al. 2003). 

4. Discussion and Conclusions 

Gil. 11-0. 12 is the first molecular filament observed to have a radial density profile 
steeper than r~ 3 at large radii r ^> r . A non-magnetic isothermal cylinder model (Ostriker 
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1964) provides a reasonable fit to both the profile and the measured constraints: submillime- 
ter flux, temperature, velocity dispersion, and core radius. The data might also be consistent 
with the generalized magnetic Stodolkiewicz model (Stodolkiewicz 1963; Nakamura et al. 
1993), whose density profile is identical to that of the Ostriker model, except that the core 
radius definition takes into account the effects of a helical magnetic field. The data are not 
consistent with the toroidally dominated regime of the FP model, which results in much 
more shallow density profiles (Fiege & Pudritz 2000). Work is in progress to provide a de- 
tailed parameter space exploration and numerical fits to the Ostriker model, the generalized 
Stodolkiewicz model and the full range of the FP model, including the steep profiled regime 
that is dominated by the poloidal field component. 

The IRDCs are thought to be a cold, dense, and quiescent class of molecular clouds 
(Egan et al. 1998; Carey et al. 1998, 2000). The dust properties measured for Gil. 11- 
0.12 between the mid-infrared and the submillimeter are consistent with this interpretation; 
however, the theoretical models of Ossenkopf & Henning (1994) require approximately 
10 5 yrs for grain coagulation to become important and thus place a lower limit on the lifetime 
of the cloud. The dynamical time of a self-gravitating non-magnetic isothermal filament, 
td ~ (Gpc) -1 / 2 , on which the filament might be expected to fragment and form stars, is similar 
to the crossing time t c ~ 2r /a, which can be calculated directly from the observations; 
t c ~ 2 x 10 5 yrs. Thus, the required coagulation time does not pose severe constraints on the 
equilibrium of the filament. It is harder to reconcile the filament's moderate central density 
n c = 1 — 4 x 10 4 cm~ 3 , which is well below the value of 10 6 cm~ 3 assumed by Ossenkopf 
and Henning in their coagulation calculations. It is possible that the dense gas in the core 
of the filament is clumped. Further studies of the correlation between 8/im extinction and 
submillimeter emission in this and other IRDCs should prove beneficial. 
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Fig. 1. — MSX 8/xm image showing the mid- infrared emission in the direction of Gil. 11. 
Note the narrow filamentary extinction lane. 
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Fig. 2. — 850/mi emission in the direction of Gil. 11. The intensity range is from -0.1 to 
0.7 Jy/beam. Note the strong correlation between the submillimeter dust emission and the 
mid- infrared extinction in Fig. 1. 
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Fig. 3. — Measured 8/xm flux versus 850/im flux at each position within the Gil. 11 filament. 
The solid line indicates the median 8/xm flux as a function of 850/im flux while the dashed- 
dotted lines give the range containing 50 percent of the 8/xm measurements. The overlaid 
dashed line represents the best-fit model (see text) with the dotted line marking the required 
8/im foreground flux /s(fg) = 2.5 x 10 -6 Wm~ 2 sr _1 . Note that the fit becomes poor at large 
850/xm flux due to the presence of a warm source at PI, producing internal 8/xm emission. 
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Fig. 4. — 850/xm emission profile across the Gil. 11 filament measured at two locations: 
(a) southern segment, (b) northern segment. In each panel, the two solid lines show inde- 
pendently the left and right profile from the ridge center. The overlaid curved dashed line 
represents the column density profile of an isothermal cylinder with r = O.lOpc (southern 
segment) and r = 0.12 pc (northern segment) after convolution by the 14.5" beam of the 
JCMT. The dash-dotted lines denote where the flux is half the peak value and the radius 
corresponding to the 14.5" beam. The dotted line in panel (a) shows a power-law relation 
with / oc r -3 , equivalent to the asymptotic isothermal profile. In panel (b), the flattening of 
one profile is due to the proximity of a second filamentary spur. 



